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Abstract 

In this paper we analyze a method of to approximation for the 
weak solutions of the incompressible magnetohydrodynamic equations 
(MHD) in unbounded domains. In particular we describe an hyper- 
bolic version of the so called artificial compressibility method adapted 
to the MHD system. By exploiting the wave equation structure of the 
approximating system we achieve the convergence of the approximat- 
ing sequences by means of dispersive estimate of Strichartz type. We 
prove that the soleinoidal component of the approximating velocity 
and magnetic fields is relatively compact and converges strongly to a 
weak solution of the MHD equation. 
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1 Introduction 



This paper is concerned with the analysis and description of a type of ap- 
proximation for the weak solutions of the 3 — D incompressible magnetohy- 
drodynamic equations (MHD). The magnetohydrodynamic system describes 
the macroscopic state of an electrically conductor, viscous and incompress- 
ible fluid and consists of the following set of equations 

^ dtu - -^Au + (u • V)u - S(B ■ V)B + Vvr + f V|Bp = 

5tB- — AB + (u-V)B- (B-V)u = 

div u = 
div B = 

endowed with the following initial data 

u{x, 0) = uo(x), B(x, 0) = Bo(x), 

where x € M^, t > 0. Here u, p and B are a-dimensional quantities cor- 
responding respectively to the velocity of the fluid, its pressure and the 
magnetic field. The a-dimensional number Rg is the Reynolds number, Rm 
is the magnetic Reynold number, M is the Hartmann number and S is de- 
fined as 5 — 



ReRm ' 

For simplicity from now on we set 

Re = Rm = S=l, p = TT + ^^. 

In particular the first equation in ([1]) is the balance of the momentum and 
the second is the equation for the induction, the equation divu = is the 
incompressibility constraint for the fluid and the equation div B = derives 
from Maxwell equations, see [20] . 

There have been extensive mathematical studies on the solutions to MHD 
equations ([1]). In particular, Duvaut and Lions [2] constructed a global 
weak solution and a local strong solution to the initial boundary value prob- 
lem. Properties of such solutions have been examined by Sermange and 
Temam in [22] • Furthermore, some sufficient conditions for smoothness 
were presented for the weak solution to the MHD equations in [T7], |18j . 
For completeness we recall here the notion of weak solutions for the system 

©• 

Definition 1.1. We say that a pair u, B is a weak solution to the MHD 
equations ([1]) if they satisfy ([1]) in the sense of distributions, namely 



lO JM.<i 



Vu • Vif + BjB,-9j(/9o- — UiUjdiipj — u ■ ) dxdt = [ uq ■ ipdx, 
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n( VB • — UiBjdiipj + ^iUjdiifj — u • ^ ) dxdt = f Bq • ipdx, 

for all if G C5^(E'^ x [0,r]), div tp = and 

divu = divB = in P'(M'^ x [0, T]). 
Moreover the following energy inequality holds 

{\u{x,t)f + \B{x,t)\^)dx + [ [ (|Vu(x,s)p + |VB(x,s)|2)dx(is 
' Jo Jr'^ 
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< J / (|uoP + |BoP)(ix, for all t > 0. 



Motivated by many applications there have been also considerable ef- 
forts to develop numerical approximation methods for the system ([1]). One 
of the major difficulty regards the development of numerical schemes in- 
cluding in an efficient way the incompressibility constraints. In fact the 
"div" operator can never be exactly zero, discretization errors accumulate 
at each iteration and after a significant amount of error accumulation, the 
approximating algorithm breaks down. In the case of the incompressible 
Navier Stokes equations in a bounded domain, Chorin [2], [3], Temam [30] , 
[3T] and Oskolkov [23], to overcome the computational difficulties connected 
with the incompressibility constraints introduced what they named "artifi- 
cial compressibility approximation". They considered a family of perturbed 
systems, depending on a positive parameter e, which approximate in the 
limit the Navier Stokes equation and which contain a sort of "linearized" 
compressibility condition, namely the following system 

dtu' + V/ = fiAu' - (u" • V) - ^(div u^)u= 

2 (2) 

edtp'^ + divu*^ = 0, 



where e is the artificial compressibility. The papers of Temam [30] , [31] and 
his book [32] discuss the convergence of these approximations on bounded 
domains by exploiting the classical Sobolev compactness embedding and 
they recover compactness in time by the well known J.L. Lions [21] method 
of fractional derivatives. The same system was used in [TH] and in [12jin the 
case of the Navier Stokes equations in the whole space ffi^ and in exterior 
domain and then modified in a suitable way in [^ for the Navier Stokes 
Fourier system in M^. As we can notice the constraint "divn = 0" has been 
replaced by the evolution equation 

dtp"^ = — - divu^, 

£ 

which can be seen as the linearization around a constant state of the con- 
tinuity equation in the case of a compressible fluid. Concerning the first 
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equation of the system ([2]) we can observe that compared to the equation 
of the balance of momentum in the incompressible Navier Stokes equations 
it has the extra term — l/2(div n^)ti^ which has been added as a correction 
to avoid the paradox of increasing the kinetic energy along the motion. In 
analogy of what has been done for the incompressible Navier Stokes we try 
to adapt this approach to the MHD equations ([1]) , keeping in mind the arti- 
ficial compressibility term we try to add to our system a similar term for the 
magnetic field, namely we introduce the following approximating system 

f 1 

dtu^ + V/ = ^Au^ - (u^ • V) u"^ - -(div u^)u^ + curlB'' x 

^ dtB' + V4>' = AB^ + curl(u^ x B") 
edtp"^ + div = 
edtct)" + div = 

where the introduced new variable c^*" is a scalar potential function. In 
the incompressible Navier Stokes equations, plays the role of a scalar 
potential that constrains the velocity to the space of divergence free vector 
fields, here in the formulation Q, the scalar potential, (j)'^ , plays the same 
role but, this time, with the aim to constrain the magnetic field, B^, to the 
space of divergence free fields. Hence, an analogous artificial compressibility 
term, which consists of the time derivative of the scalar potential function, 
(/>^, has been added to the solenoidal constraint equation, and similarly, a 
Vt/)'" term with some other stabilization terms are included in the induction 
equation. A similar approximating system has been introduced in [33] and 
[34| for a MHD system without viscosity but the analysis is performed only 
from a numerical point of view. In this paper we analyze rigorously the 
convergence of the system ([3D towards the system ([T]) as e — in the case 
of the whole space . Furthermore we assign to the system Q the following 
initial conditions 

u%x, 0) = ug(x), B%x, 0) = Bg(x), 
/(x,0) =pg(x), 0^(x,O) =</.g(x). 

The limiting behavior as e J, of the initial data ([H deserves a little discus- 
sion. Indeed Q requires the initial conditions Q while the MHD equations 
([!]) require only the two initial condition for the velocity u and for the 
magnetic field B. Hence, our approximation will be consistent if the ini- 
tial datum on the pressure p^ and of the potential <j)^ will be eliminated 
by an "initial layer" phenomenon. Since in the limit we have to deal with 
weak solutions in the sense of the Definition 11.11 it is reasonable to require 
the finite energy constraint to be satisfied by the approximating sequences 
{u^ ,B^ ,p^ , (j)^). So we can deduce a natural behavior to be imposed on the 
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initial data (uq, Bq,Pq, (/)q), namely 

= u^(-,0) — ^ uo = u(-,0) strongly in L^(R^) ^ 
Eg = B^(-, 0) — ^ Bo = B(-, 0) strongly in L'^{R^ 
= i/F/(-,0) — y strongly in L^(R^) 
= Ve0^(-,O) — > strongly in L^{W^). 



(ID) 



Let us remark that the convergence of ^/epQ and -v/e^g *o ^ is necessary to 
avoid the presence of concentrations of energy in the limit. Now we can 
state our main result. 

Theorem 1.2. Let (u^, B^,p^, 0'^) be a sequence of weak solution in of 
the system assume that the initial data satisfy (ID). Then 

(i) There exist u, B G L°°([0, T]; L2(r3)) n L2([0, T]; i/^ (M^)) such that 

u^^u weakly in L^{[Q,T];H^{W^)). 
B^ ^ B weakly in L^{[Q,T\,H^{R^)). 

(ii) The gradient components Qu^ , QB^ of the vector fields u^, B"^ satisfy 
respectively 

Qu^ — > strongly in L2([0, T]; LP(]R^)), for any p G [4,6). 
QB= — > strongly in L^([0, T];LP{R^)), for any p G [4, 6). 

(iii) The divergence free components Pu'^ , PW of the vector fields u^, B"^ 
satisfy 

Pu^ — > Pu = u strongly in L'^{[0,T]; Lf^^{R^)). 
PB^ — ^ PB = B strongly in L'^{[0, T]; LI^{R^)). 

(iv) (u = Pu, B = PB) is a Leray weak solution to the incompressible 
magnetohydrodynamic system 

P{dtu - Au + (u • V)u - (B • V)B) = 0, 

dtB - AB + (u • V)B - (B • V)u = 0, 

and the following energy inequality holds 

- [ (|u(x,t)p + |B(x,t)|2)(ix+ / / (|Vu(x,t)|2 + |VB(x,t)|2)dxdi 
< \ f (|u(x,0)|2 + |B(x,0)|2)dx. (5) 
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One of the major issues connected with the proof of such kind of result 
is the presence of acoustic waves for the pressure and for the potential 
(f)'^ which exhibit fast oscillations when e — )■ 0. As we can see by considering 
together the equations and ([3]) 3 or ([3]) 2 and (IS])^ as e goes to at the 
level of the pressure or of the potential, the acoustic waves propagate with 
high speed ^ in the space domain. Because of the fast propagating of the 
acoustics one expects the velocity and the magnetic field to converge 
only weakly. In fact if, for instance, we project the third equation of ([3]) on 
the space of divergence free vector fields we get that 

dtPW = P(-AB^ + curl(u^ x B=)), 

in other words the soleinodal component PB^ is relatively compact in the 
time variable. Therefore time oscillations are supported by the gradient part 
of the field B^ and are related to the propagation of acoustic waves. We 
overcome this trouble by using the dispersion of these waves at infinity yield- 
ing the strong convergence of and B^ in space and time. In particular the 
system ([3]) will be discussed as a semilinear wave type equation for the pres- 
sure function and for the potential function and the dispersive estimates will 
be carried out by using the L^^-type estimates due to Strichartz [16} [T9t [29]. 
The connection with the dispersive analysis of the acoustic wave equation 
has also been considered to study the incompressible limit problem. Similar 
phenomena appear also in the modeling the Debye screening effect for semi- 
conductor devices, [11], [13]. In particular a similar wave equation structure 
has been exploited in various way by the paper of P.L. Lions and Masmoudi 
[22], Desjardin, Grenier, Lions, Masmoudi [6], Desjardin Grenier [5] -It is 
worth to mention here that this type of singular limits from hyperbolic to 
parabolic systems is not covered and doesn't fits in the general framework 
of diffusive limits analyzed in [9j. 

This paper is organized as follows. In Section 2 we recall the notations, 
some mathematical tools needed in the paper and recall same basic defini- 
tions. In Section 3 we recover all the uniform bounds that can be deduced 
by the only use of the energy inequality. Section 4 is devoted to the study 
of the acoustic equations. In Section 5 we prove the strong convergence of 
the velocity and magnetic field. Finally in Section 6 we give the proof of the 
main result. The last Section 7 is devoted to the study of the artificial com- 
pressibility approximation for general unbounded domain, as for example an 
exterior domain. 

2 Notations 

For convenience of the reader we establish some notations and recall some 
basic facts that will be useful in the sequel. 
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We will denote by V{R'^ x M+) the space of test function C^(R'^ x M+), 
by P'(M'^ X the space of Schwartz distributions and (•,•) the duality 
bracket between V and V and by MtX' the space C°([0, T]; X)'. Moreover 
W'''P{R'^) = (/ - A)-tLP(]R'^) and = W^'='2(M^) denote the non- 

homogeneous Sobolev spaces for any 1 < p < oo and A; G M. VF'^'P(R'^) = 
(I- A)- ^2 LP {M.'^) and = M^^'2(R'^) denote the homogeneous Sobolev 

spaces. The notations L^Ll and L^Wx''^ will abbreviate respectively the 
spaces LP([0,r];L9(M°')), and U'{[{),T]]W^^''{W^)). 

We shall denote by Q and P respectively the Leray's projectors Q on the 
space of gradients vector fields and P on the space of divergence - free vector 
fields. Namely 

Q = VA"^ div P = I-Q. (6) 

Let us remark that Q and P can be expressed in terms of Riesz multipliers, 
therefore they are bounded linear operators on every W'^'P (1 < p < oo) 
space (see [28]). 

We recall the following vector identities that we will use later on in the 
paper. If A and H are two vector fields then it holds 

V(|A|2) = 2(A- V)A + 2A X curlA, (7) 

curl( A X H) = A div H - H div A + (H • V) A - (A • V)H, (8) 
div((A X H) X H) = curlH x H • A - curl(A x H) • H. (9) 

3 Uniform bounds 

In this section we wish to establish the priori estimates independent on e, 
that comes directly by the energy associated to the system ([3]). In fact we 
define the energy of the system ^ as 

m = }; [ (|u^(x,t)P + |B^(x,t)|2 + |VF/(x,t)P + |V^0^(x,t)|2)dx. 

(10) 

If we multiply the equations of the system ([3]) respectively by u^, B^, and 
(j)^, we sum up and integrate by parts in space and time, with the use of the 
identities ([T])-® we are able to prove the following theorem 

Theorem 3.1. Let us consider the solution (u^, B^,p^, c/)^) of the Cauchy 
problem for the system Assume that the hypotheses (ID) hold, then one 
has ^ 

^(*)+ / / (|Vu^(rE,s)p + |VB^(x,s)p) dxds = S(0). (11) 

Jo JR3 
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(15) 



Prom the energy equality (jlip we get the following list of estimates 

Ve/, ^fe(\f are bounded in L°°([0, T]; l2(m3)), (12) 
Vu^ VB^ is bounded in ^^([0, T\ x M^). (13) 
By combining ()12p with Sobolev embeddings theorems we obtain that 

u^B^ are bounded in L°^([0,T];L2(M3))nL2([0,T];L6(M3)), (14) 

Using together (fT3]) and ([H]) we have that 

(u^-V)u^,u^divu^,u^ X B^,curlB^ x B^ 
ar e b ounded in ( [0 , T] ; L ^ (M^ ) ) n L 1 ( [0 , T] ; L^/^ (^3 ) ) ^ 

Finally combined with the equations ([3]) 3 and ([3]) 4 gives rise to 

epl,£(j)l relatively compact in i7"^([0,T] X M^). (16) 

4 Acoustic wave equations and estimates 

As already mentioned in the Introduction one of the major issue in this 
type of analysis is the presence of acoustic waves which prevents the strong 
convergence of the gradient part of the velocity and magnetic field. It turns 
out that these waves satisfy inhomogenous wave equations which can be 
seen as the superposition of waves with different frequencies scales connected 
with the viscous tensor of the fluid and to the convective terms. In order to 
study this phenomena first of all we rescale the time variable, the velocity, 
the magnetic field, the potential 0^ and the pressure in the following way 

t 

iL{x,T) = U^{x, ^/et), B{x,t) ='B^{x,^/£T) 

p{x,t) = p''{x,y/eT), (t>{x,T) = (|)%X,^/eT). 



(17) 



As a consequence of this scaling the system ([3]) becomes 

/s 

9t-u + -y/eVp = -y/eAu — -y/e (u • V) u 2"^*^^^ ^ B X B 

drB + ^/eV4> = ^/eAB - ^curl(ii x B) 
^/edrP + div u = 
Ve9^(^ + divB = 

(18) 

then, by differentiating with respect to time the equation (jlSP o and by using 
(jlSP ]^. we get that p satisfies the following wave equation 



drrP - Ap= -Adivu + div ( (u • V) u+ ^(divu)u + curlB x B j , (19) 
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where, by using the bounds ()12p - ()15p . the inhomogenous parts belongs to 
L2([0, T];H-^{M.^)) + Li([0,T]; I^-1'3/2(R3)). 

In the same way by differentiating with respect to time the equation ()18|) ^ 
and by using (fT8|) n. we get that B satisfies the fohowing wave equation 

- = -AdivB, (20) 

with the right-hand side in L^([0, T]; /^"^(M^)). In order to estimate p and 
(j) we will make use of the dispersive estimates of Strichartz type that we 
briefly recall in the next section. 

4.1 Strichartz estimates 

If u; is a (weak) solution of the following wave equation in the space [0, T] xM'^ 

g + A) w{t,x) = F{t,x) 
[wiO,-) = f, dM0,-)=9, 

for some data f,g,F and time < T < oo, then w satisfies the following 
Strichartz estimates, (see [TB], [TP] ) 

IHlL^Lr + \\dtW\\ ^,^-l,r < WfWjjy + + \\F\\^f^,J, (21) 

where {q,r), {q,r) are wave admissible pairs, namely they satisfy 

2 /I 1\ 2 /I 1 

q-^ ~ '\2~r) '^-^ ~ \2~7 

and moreover the following conditions holds 

1 d d 1 d „ 

q r 2 q r 

In particular, later on we shall use (j2ip in the case ofd = 3, (g', f') = (1,3/2), 
then 7 = 1/2 and {q,r) = (4,4), namely the following estimate 

MlI^ + \\dtw\\^,^-iA < 11/11^1/2 + + WHlilI^^- (22) 

Beside the Strichartz estimate ()2ip or ()22p in the case of d = 3 (see [27j). 
by using Duhamel's principle and an earlier Strichartz estimate |29] , we can 
also deduce the estimate 

\\w\\l4^ + \\dtw\\ ^^^^^1,4 < 11/11^1/2 + ll5ll^-i/2 + II-^IIlili- (23) 
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4.2 Estimates for (f> 

In order to apply the dispersive estimates of the type ()2ip to the wave 
equation (|20p in different components we have to specify the initial data, 
namely we will analyze the following wave equations 

(^(x, 0) = (9^(/)(x,0) = divB^, 
We are able to prove the following theorem. 

Theorem 4.1. Let us consider the solution (u*", B^,p^, </)^) of the Cauchy 
problem for the system Assume that the hypotheses (ID) hold. Then we 
set the following estimate 

+ \/r||divB^||^2i2. (25) 



Proof. Since (f) and is a solution of the wave equation (|24|) and G is bounded 
in L\H~'^ we can apply the Strichartz estimate ()23p . with {x,t) G R'^ x 
(0,T/y^) and we get 

I|A-^0||l4^^ + ||5,A-l,^||i4H.-l,4 < ||A-10(X, 0)11^1/2 + ||A-V2a^^(^, 0)11^-1/2 

+ ||A-iG||^,^3/2 (26) 



namely 



+ W^AWl^^jv-^:^ ^ [[-^(a;, 0)11^-3/2 + ||(9r<^(x, 0)11^-5/2 



-2,4 



+ -^||divB||^2^2 (27) 



L'-((0,T/v^);L9(I83)) — £ ||(^ IIL'-([0,T];L9 



Finally, since 



we end up with (j25|) . □ 
4.3 Estimates for p 

In order to get the estimates for p we split it as p = pi + p2 where p\ and 
P2 solve the following wave equations: 

[drrVx - Api = -A div u = Fi ^^^^ 
|pi(rE,0) = (9^pi(x,0) = 0, 
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9rrP2 — Ap2 = div ^(u • V) u + ^(divu)u + curlB x = F2 ^^^^ 

^P2{X,Q) =p{x,0) drP2{x,0) = drPix,0). 

Now by following the same line of arguments as in the previous section or 
as in [To] we can prove the following theorem. 

Theorem 4.2. Let us consider the solution (u^, B'^,p'^, (^^) of the Cauchy 
problem for the system Assume that the hypotheses (ID) hold. Then we 
set the following estimate 

^'^'lIP'llLfH/r^.* + e'/'l|9iP"lli4^-3.4 < V^MWlI + II div ll^-i +^/T|| div 11^2^2 

+ II (u^ • V)u^ + ^(divu^)u^ + curlB^ x B^|Li^3/2. 

(30) 



5 Strong convergence 

This section is devoted to the strong convergence of the magnetic field B^ 
and the velocity field u^. In particular we will show that the gradient part 
of the magnetic field QB^ and of the velocity Qu^ converges strongly to 
0, while the incompressible component of the magnetic field PB^ of the 
velocity field Pu^ converges strongly to PB = B and Pu = u respectively , 
where B is the limit profile as e 4, of B^ and u is the limit profile as e 4- 
of u^. We start this section with a short remark on some strong convergence 
which are an easy consequence of the estimates (fT^ . ([2^ and (pin|) . whose 
proof is left to the reader. 

Proposition 5.1. Let us consider the solution (u^, B^,p^, (/)^) of the Cauchy 
problem for the system ([3]). Assume that the hypotheses (ID) hold. Then, 
as e 10, one has 

strongly in L°^{[0,T];L^{M.^)) D L^{[0,T];W~^^^iM.^)), 
strongly in L°°{[0,T]; L'^{R^)) n L^{[0,T];W~'^'^{R^)), 
strongly in W^^''=^i[0,T]; L'^{M.^)) D L^{[0,T];W^^'^iM:^)). 
strongly in W-^'^{[0,T]; L^ (M^)) n L"^ ( [0, T] ; W'^^^ (M^)) . 

Before going into the analysis of these strong convergences we list here 
some technical results which will use later on. 

Lemma 5.2. Let us consider a smoothing kernel ip G Cq°{M.'^), such that 
■0^0? /]Rd i^dx = 1, and define 



ep" — > 
e(j)'' — > 
div — y 
div B^ — >0 
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Then for any f € H^{W^), one has 

ll/-/*V'a||LP(M.) <CX""I|V/||l2(k.), (31) 

where 

pG[2,oo) ifd = 2, ps[2,6] if d = 3 and a = d(- 



2 p 



Moreover the following Young type inequality hold 



\\f*i'a\\LP(R'i)<Ca " '^Vi (32) 



for any p,q € [l,oo], q < p, s > 0, a £ (0, 1). 

Proposition 5.3. Let he T <Z L^([0, T]; 5), 1 < p < oo, B a Banach space. 
J- is relatively compact in L^([0,T];i?) for 1 < p < oo, or in C{[0,T]; B) 
for p = oo if and only if 



(i) 



J f{t)dt, f €z B^ is relatively compact in B , < ti < t2 < T, 



(ii) lim \\f{x + h) - f{x)\\LpnQT-h\,B) = uniformly for any / G J". 



5.1 Strong convergence of QW and FB^ 

Proposition 5.4. Let us consider the solution (u"^, B^,/>^, (/)^) of the Cauchy 
problem for the system ([3]). Assume that the hypotheses (ID) hold. Then as 

— > strongly in L'^{[Q,T]]L^{W^)) for any p £ [4,6) . (33) 
Proof. In order to prove the convergence ([33|) we split QB"^ as follows 

\\QB'\\l2l. < WQB' - QB' * Valli^Lg + IIQB^ * ^alL^Lg = Jl + J2, 

where ipa is the smoothing kernel defined in Lemma |5.2[ Now we estimate 
separately Ji and J2. For Ji by using ([3T]) we get 



Ji<a'"V2 p; ||VQB-'(i)|||2dij < a' "V2 ||VB"||^2i2 . (34) 

Hence from the identity QB^ = —e^^^V /S.^^ dtcjf and by the inequality 
([32]) we get J2 satisfies the following estimate 



(35) 
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Therefore, summing up (|34p and ()35p and by using ()13p and (|30p . we con- 
clude for any p G [4, 6) that 

Finally we choose a in terms of e in order that the two terms in the right 
hand side of the previous inequality have the same order, namely 

a = e^/is. (37) 

Therefore we obtain 

\\QB'\\l2lp < Cte^ for any p e [4,6). 

□ 

It remains to prove the strong compactness of the incompressible com- 
ponent of the magnetic field. To achieve this goal we need to look at some 
time regularity properties of PB^ and to use the Proposition 15.31 

Proposition 5.5. Let us consider the solution (u^, B'^,p'^, c/;^) of the Cauchy 
problem for the system ([3]). Assume that the hypotheses (ID) hold. Then as 

PB^ — > PB, strongly in ^^(0, T; Lf^^(R^)). (38) 

Proof. From (|13p an (|14p we know that PB^ is uniformly bounded in L/^H].. 
The strong convergence ([38|) follows by applying the Proposition 15.3^ pro- 
vided that for all h G (0, 1), we have 

||PB^(t + /i)-PB^(t)||^2([o,T]xR3) <C7t/i^/^ (39) 

The rest of the proof is devoted to show (f39|) . Let us set = B'^{t + h) — 
B^[t), we have 

|PB^'(t + h)- PB^(t)||i.(p^^],j,3) = j'^ j^^ dtdxiPZ') ■ (PZ^ - PZ^ * V'a) 

+ / / dtdx{PZ') ■ (PZ^ * Va) = h+ h- 

Jo Jw^ 

(40) 

By using ()3ip we can estimate Ii in the following way 

Ji < WPZ'U^Ll r \\PZ%t) - iPZ' * ^|Ja)mLldt 
Jo 

<«rV2||B-||,c„^.||VB^||^2 . (41) 
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Let us reformulate PTf in integral form by using the equation ([3]) 2, hence 



li < 



V2' 

rT f ft+h 

dt dx ds[^W - curl(u^ x B^))(s, x) ■ [PTI" * ijc,){t, x) 
Jo Jt 



■ (42) 



Then integrating by parts and by using (|32p . with p = 00 and q = 2, we 
deduce 

/ rt+h \ 1/2 

h < h\\VB^\\l2 + Ca~^/'^T^/'^\\B^\\Lo.^2 ih / ||u" x B^Wl^ds \ 

< /i||VB^||2,^ +Ca-5/2rV2/,||B-||^c«i2||u- X B^W^^^.. (43) 

Summing up Ii , I2 and by taking into account (fT3]) , , (fT5]) , we have 

and, by choosing a = h?/"^ , we end up with ()5U]) . □ 
5.2 Strong convergence of Qu^ and Pu^ 

The strong convergence for Qu^ and Pu^ can be performed with the same 
method as in the previous section. In particular we are able to prove the 
following propositions. 

Proposition 5.6. Let us consider the solution (u^, B'^,p'^, (/)^) of the Cauchy 
problem for the system ([3]). Assume that the hypotheses (ID) hold. Then as 
eiO, 

Qu^ — > strongly in L'^{[0, T];LP{R^)) for any p £ [4, 6) . (45) 



Proposition 5.7. Let us consider the solution {u'^ ,B^ ,p^ of the Cauchy 
problem for the system ([3]). Assume that the hypotheses (ID) hold. Then as 

Pu^ — > Pu, strongly in L^{0, T; L^„^(M^)). (46) 
For more details on the proof of the Propositions 15.61 and 15.71 see also 

m, m- 



6 Proof of the Theorem 11.21 

The weak convergence in (i) is a consequence of the uniform bounds (|13p 
and (|14p . The strong convergences in (ii) and (iii) follows from Propositions 
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\5.7\ 15.51 \5.6\ 15.41 In order to prove (iv) we apply the Leray projector P to 
the equations (l3|)^ 21 then, taking into account the identities ©-([SI) it follows 

atPu^-/iAPu^+P((u^- V)u=) + ^P((divu^)u^)-P((B^-V)B=) = 0, (47) 

dtPB^ - APB" + curl(u^ x B^) = 0. (48) 

Now, by testing the equations (|47|) and (|48|) against a test function it is easy 
to pass into the limit in the linear terms. For what concerns the nonlinear 
terms it is enough to decompose the vector fields and B^ in their diver- 
gence free part and gradient part end then to pass into the limit by using 
(ii) and (ill). The nonlinear term curl(u^ x B^) in the equation ()48p can 
be handled by combining the previous decomposition with the identity ([8]). 
For more details in this step see [10], [12], [13] . 

Finally we prove the energy inequality ([5]). By using the weak lower 
semicontinuity of the weak limits, the hypotheses (ID) and by denoting by 
X the weak-limit of y/ep'^ and by S the weak-limit of -v/e</>^ , we have 

I [ {\x\' + \^\')dx 
+1 [ (|u(x,t)|2 + |B(x,t)p)dx+ / / (|Vu(x,t)|2 + |VB(x,t)p)(ix(it 
<liminffi/ ^\^p^\^ + \^^^^)dx + l [ i\u^{x,t)\'' + \B^ix,t)\^)dx) 

+ liminf / [ {\Vu'{x,t)\'^ + \VB'{x,t)\'^)dxdt 
= liminfi/ (|u§|2 + |B^ol' + I^^Pol' + l^^'ASl')d^ = ^ / (|uoP+|Bo|2)dx 

e^O 2 7iR3 ' '2 J^3 

7 General unbounded domain case 

This last section is devoted the analysis of the artificial compressibility ap- 
proximation in the case of a general unbounded domain as for example the 
half-space, an exterior domain, an unbounded strip, etc. For what concerns 
the proof of the energy estimate (|lip and of the strong convergence of the 
soleinoidal components Pu"^ and PB^ of the velocity and magnetic field we 
can easily adapt the proof in this case. The main difference is in the method 
used to prove the compactness of the gradient component of and B^ . As 
already pointed out the main issue is to control the time oscillations due to 
the presence of acoustic waves. Since we are in an unbounded domain also 
in this case we have a dispersive phenomena for the acoustic waves but it 
cannot be handled in the same way as in the previous part. In fact, as is 
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well known, the Strichartz estimates in unbounded domain different than 
M*^ become more delicate and usually require some restrictions on the shape 
and geometry of the boundary. For instance in the papers [1], [23], [26] the 
domain is required to be star-shaped or non trapping. Here we carry out an 
approach which can be applied to any unbounded domain il. which satisfy 
the hypotheses 

• the point spectrum of the Neumann Laplacian A^r in L^(r2) is empty. 

Although the absence of eigenvalues for the Neumann Laplacian represents, 
in general, a delicate and highly unstable problem there are numerous exam- 
ples of such domains - the whole space M^, the half-space, exterior domains, 
unbounded strips, tube-like domains and waveguides. Similar approach has 
been used also in [15], [8]. 

Here, for simplicity we will perform our analysis for an exterior domain 
(the proof can be easily adapted to other unbounded domains) end we will 
show the convergence for the gradient part of the velocity field u^. The 
convergence for can be handled in the same way. 

We say that is an exterior domain if it is the complement in of a 
compact set (usually called compact obstacle). Since Q, has a boundary dil, 
we have to be careful in defining the Leray projectors. In particular for a 
vector field v we define the Helmholtz decomposition as 



divPv = 0, Qv = Vf, Alp = divv, • n\gn = v • n\dn, 

and n is the exterior normal to the boundary 9(7. We consider the sys- 
tem ([3]) endowed with the initial conditions ppp . we specify the following 
asymptotic behavior of and B"^ 



V = Pv + Qv 



(49) 



where 



u'' ^ B^ ^ 



as Ixl — )■ 



We assume that the boundary is impermeable, namely 



• n\QQ = 



B^ 



nbn = 



and the no-slip boundary conditions hold 



B^ 



an = 0. 



First of all we rewrite the equations and (I3|)3 in terms of the acoustic 
potential 'f^. 




(50) 
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where F*^ = /iAu^ - (u- • V) - ^ (div u^)u= + curl x and V/= = QF^ 
The sohitions for ()50p may be expressed by means of Duhamels formula as 



^%t) = ^exp j^iV-ATv- 



+ 2 exp I -IV-Aat- 



^^(0) + 



i 



/(O) 
/(O) 



(51) 



N 



+ 



-Aiv^— ^ ) + exp ( -i-v/-AAf 



t - s 



[ns)]ds, 



where A^r is the Laplace operator endowed with the homogenous Neumann 
boundary conditions in the Hilbert space L^(r2). By observing the formula 
()5ip we see that oscillates fast in time at the frequency 1/e as soon as 
— Ajv possesses a positive eigenvalue. Since we are in an unbounded domain 
— Ajv has an empty point spectrum, this allows us to apply the celebrated 
RAGE theorem, see Cycon et al. [3](Theorem 5.8), see also [8]: 

Theorem 7.1. Let H be a Hilbert space, A : 'D{A) C H ^ H a self-adjoint 
operator, C . H ^ H a compact operator, and Pc the orthogonal projection 
onto the space of continuity He of A, specifically, 



H = Hr 



clH|span{ 



w & H \ w an 



eigenvector of A}^. 



Then 



exp(-itA)CPcexp(it^) dt 



as T —)• oo. 



CiH) 



(52) 



We apply Theorem ITU to H = L'^{n), A = ^/^^A^, C = x^G{-An), 
with X e C^{^), X > 0, G e C^(0,oo), < G < 1. Taking r = 1/e in 
we have 

£ ^exp (^-iiy^^^ x'G(y^A^)exp j^i^T^A^^ X;y'^ dt 

< ^(^)ll^llL2(n)ll^llL2(n), 
where uj{e) — )• as e — )■ 0. If we take Y = G{\/—A]\f)[X] we deduce that 

r xG(y^)exp (il^/^) [X] dt < a;(e)||X||2 (53) 

for any X £ L^(f]), for more details see [8]. The relation (j53p implies the 
local decay to zero of namely 

I|V^'^IIl2(E') for any compact set K C (0,r) x $7, 

hence the strong compactness in space and time for the gradient part of u^. 
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